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ABSTRACT 

We quantify the importance of the mechanical energy released by radio-galaxies inside galaxy groups. 
We use scaling relations to estimate the mechanical energy released by 16 radio- AGN located inside 
X-ray detected galaxy groups in the COSMOS field. By comparing this energy output to the host 
groups' gravitational binding energy, we find that radio galaxies produce sufficient energy to unbind 
a significant fraction of the intra-group medium. This unbinding effect is negligible in massive galaxy 
clusters with deeper potential wells. Our results correctly reproduce the breaking of self-similarity 
observed in the scaling relation between entropy and temperature for galaxy groups. 

Subject headings: galaxies: clusters: general — intergalactic medium — radio continuum: galaxies — 
X-rays: galaxies: clusters — galaxies: active 



1. INTRODUCTION 

Galaxy groups are important laboratories in which 
to investigate the importance of non-gravitational pro- 
cesses in structure formation. These processes are poten- 
tially more important in galaxy groups than in massive 
clusters because of their lower gravitational binding en- 
ergy. This is suggested by the significant deviation of 
the observed X-ray luminosity and entropy versus tem- 
perature {Lx-T and S-T) scaling relations in groups 
compared to the relatio n expected in a purely gravita- 



tional scenari o (see also [P ratt fc Arnaud'^ 2QQ3 



Ma rke- 

vitch 1998; Arnaud & Evrard 1999; Ponman et dIF' 



Sun et al.||2QQ9[ [Pratt 'eral.||2QQ9D . Radiative cooling 
can be invoked to explain this deviation, but then the 
predicted fraction of stars in clusters of a given mass is 
incorrect ( |Voit||2QQ5[ [Balogh et al ||2QQ8D . To simultane- 
ously explam the properties of the mtra-cluster /group 
medium (ICM) and account for the observed properties 
of galaxies, it is necessary to take into account a major 
contribution to the cluster /group energetics from non- 
gravitational heating. 
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The two main sources of non-gravitational heating are 
star-formation and active galactic nuclei (AGN). Cos- 
mological simulations (e.g. Kay||2QQ4 Bower et al.||2QQ6 
^Sijac ki & Springel 20061) snow that both processes are 
required to reproduce the properties of the ICM. In par- 
ticular, recent simulations by Bower et al.| (2008) suc- 
cess fully reproduce both the galaxy and ICM properties 
(see Short & Thomasj2009j) when they include a "radio- 
mode" AGN feedback phase: in this phase the movement 
of bubbles inflated by the AGN jets transfers energy into 
the gas within the cluster (mechanical heating). The ob- 
servable objects providing this type of feedba ck inside 
groups and clusters would be radio galaxies (iCroton et 
al. 2006). The main difference between the Bower et al. 
( 2008 ) model and o thers, including radio-mode AGN 



(.Bower et al. 2006 Sijacki & Springel 



2006 



et al.||2008[ ) is that it allows the radio mode 



Puchwein 
feedback to 

expel gas from the X-ray emitting regions of the system. 
The importance of such AGN-fee dback in g roups could 



Mc- 



explain the observational result bylLin et al. ( 2003p 
,Carthy et al., (2007) and Giodini et al.| (200g)7that^tlie 
total baryon fraction in groups is lower than the cos- 
mic value estimated fro m cosmic microwave background 
(CMB) observations (see [Giodini et al.|2009|for more de- 
tails). The discrepancy decreases in systems of higher 
total mass, such that it is <l<j for massive clusters. 

In this Paper we propose a simple, direct method to 
test the hypothesis that radio galaxies in groups can 
indeed inject enough mechanical energy to unbind the 
intra-cluster gas. The Paper is structured as follows. In 
Section 2 we select a sample of 16 groups from the COS - 



MOS 2 deg^ survey discussed in [Giodini et al. 
each hosting a radio galaxy within the viria. 



(2009|, 
radius 



eacn nostmg a radio galaxy witmn tne virial radim 
dSchinnerer et al.[|2007[ [Smolcic et al 2008) , plus a con- 
trol sample of massive clusters from Birzan et al. (2004) 
In Sections [3] and [5] we then compare the groups bind- 
ing energy to the mechanical energy output by the ra- 
dio sources, derived from their total radio luminosity 
through scaling relations. Applying this method, we 
show that the mechanical removal of gas from the group 
region is indeed energetically feasible for systems below 
^SxIO^^Mq. In Section [6] we discuss how this scenario 



2 



Giodini et al. 



compares to the deviation in the scahng relation between 
entropy and temperature at the groups scale. 

We adopt a ACDM cosmology with /i=0.72, 1^^=0.25, 
1]a=0.75. 

2. THE SAMPLES 

2.1. Radio galaxies in X-ray detected groups 

We use the catalog of 91 X-ray se lected groups from 
the COSMOS survey ( [Scoville et aL| [2QQ7a ; Finoguenov 
et al. in p reparation), selected as described in Giodini 
et al. (2009j). Extended source detection was performed 
usmg a multiscale wavelet reconstruction of a mosaic 
of XMM and Chandra data. For each group, member 
galaxies are identified within Rsoo^^ of the group center, 
utilizing the high quality photometr ic redshifts availab le 



Ilbert et al. 



'Smolcic et al 



2009). 



(cr(Az)/(l + z)=Om at lAB^^ <25, 
We use a sub-samp l e of t he V 
log (Schinnerer et al. [2007 
identify radio galaxies lying inside the X-ray selected 
groups. Of the 60 radio galaxie s identified within the 



cata- 

2008^ 



to 



VLA-COSMOS Large Project ( [Schinnerer et al.||2007 
1.49 GHz), about 80% have been associat ed with a secure 
optical counterpart (Smolcic e t al.|[2Q08 ) with zab ^ 26, 
and accurate phot ometry (thus also with accurate pho- 
tometric redshifts; [Ilbert et al |2009[ [Salvato et aT|2 009 ) . 

We have cross-correlated this sample of radio galax- 
ies with the X-ray selected galaxy groups in 3D space 
using a search radius of 1 x R200 (Finoguenov et al. 
in preparation) around the groups' centers and within 
0.02x(l+z) from the group's redshift. This resulted in 
a sample of 16 systems matched in position and red- 
shift. In Appendix A we show the contours of the 
radio 20 cm and X-ray emission superimposed to the 
SUBARU zp band image for each of the groups. In 9 
out of 16 cases the radio galaxy is located in the core 
of the group (defined as R <0.15i?20o)- The 20 cm ra- 
dio luminosity densities of these galaxies range from 
^ 5.5x 10^^-4.8x 10^^ W Hz"^, with a median luminos- 
ity of 8.9x10^"^ WHz-i(i/F^ - 7.3x10^^-6.4x10"^^ 
s~^, with a median luminosity of 1.18 x 10^^ erg s 
. This median luminosity is at the high end of the ra- 
dio luminosity distr ibution of the full ra dio AGN sam- 
ple (c.f. Fig. 17 in [Smolcic et al. 2008 and Fig. 5 in 



erg 
-1^ 



Smolcic et al. 



j2009[ ), consistent with previous findings 
that powerf ul radio galaxies in habit group- scale environ- 
ments (e.g. [ Baum et al.[[1992 ). The redshift distribution 



of the 16 groups is fairly uniform between 0.1 and 1, with 
the exception of 6 sources concentrated at z~0.3 (where 
a large structure extends throughout the whole COS- 
MOS field). The groups have X-ray luminosities ranging 
from 1x10^^ to 8.7x10^^ erg s~^ and span a mass range 
of 2x10^^ <M2oo <2xlO^'^M0 with a median mass of 

Ra (A=500,200) is the radius within which the mass density 
of a group/cluster is equal to A times the critical density (pc) of 
the Universe. Correspondingly, Ma = A pc{z) {4:7r/3)R^ is the 
mass inside Ra- M200 is computed using an LX-M200 relation 
established via the weak lensing analysis in Leauthaud et al. 2010 
The catalogue value of M200 is converted into M500 assuming an 
NFW profile with a concentration parameter computed from the 
mass-dependent relation of Maccio et al. (2007). 

AB magnitude in the SUBARU i band. 

The term "radio galaxy" is used here to describe an extended 
radio source with clear jet/lobe structure. 

Computed using the total flux densities (Fi^). K-correction 
is also applied assuming a spectral index of a = 0.7 (Fjy oc h'~^). 



7.14xlO^^M0. 

2.2. The comparison sample of massive clusters 

The COSMOS X-ray sample is mostly composed 
of groups. We complement it with 12 well known 
radio galaxies inside massive clust ers, extracted from 
the sample of [Birzan et al. ( 2004 ) . We use those 
clusters from the Birzan 's sample which overlap with 
the HIFLUGCS survey ( [Reiprich fc Boehringerl[2QQ2 ) 
so that we can use the X-ray parameters determined 
from the HIFLUGCS clusters. In addition we require 
that the radio source within those clusters is associated 
with a secure NIR cou nterpart in the 2MASS catalogue 
(Skrutskie et al. 2006). These requirements eliminate 4 
of the clusters in the original Birzan et al. sample. Each 
of these clusters contain X-ray cavities associated with 
radio bubbles likely connected with AGN activity of the 
central galaxy. The radio galaxies have been identified 
within the NRAO V LA Sky Survey (NVSS) at 1.49 GHz 
([Condon et al.'1998 ), except in the case of the Centaurus 
Cluster where data come from the 1.41 GHz Parkes 



Radio Sources Catalogue ( [Wright fc Otrupcek 1990). 
The 20 cm radio luminosities of the radio galaxies range 
between 2x10^^ and 2 xlO^^ erg s~^, with a median 
luminosity of 1.4x10^^ erg s~^, more than 10 times 
higher than the median radio luminosity of the radio 
galaxies in the COSMOS sample. 

The X-ray parameters for these clusters are provided 
by the X-ray analysis in the HIFLUGCS survey, and 
converted for the standard cosmology used in this paper. 
The sample consists of very local clusters, ranging in 
total mass between 1x10^"^ <M2oo <1.2xlO^^M0 with 
a median mass of 4.25xlO^^M0, almost 10 times 
higher than the median total mass of the systems in the 
COSMOS sample 

3. ANALYSIS OF COSMOS GROUP SAMPLE 

3.1. Mechanical energy input by radio galaxies in groups 

We estimate the mechanical energy input by a radio 
galaxy into the ICM over the group lifetime from the 
mechanical luminosity of the radio source multiplied by 
the fraction of time a massive galaxy spends in the radio- 
AGN phase. The mechanical luminosity for the radio 
galaxies in our sample is estimated from the scaling re- 
lation presented in Birzan et al. (2008) . These authors 
studied a sample of galaxy clusters showing signatures of 
cavities and bubbles in the X-ray surface brightness 2D 
distribution, with a powerful radio source as a central 
galaxy. The cavity power of the radio source, estimated 
from the pdV work of the jet/lobe on the surrounding 
ICM, is found to be correlated (albeit with a large scat- 
ter) with the monochromatic radio power at 1.49 GHz of 
the central galaxies {Pia9GHz) as 

p p0.35±0.07 /-jx 



(see Eq. 16 in Birzan et al. 2008). Piaoghz is computed 
from the radio-emission of the entire source. This esti- 
mate is a lower limit to the mechanical luminosity of the 
AGN outbursts, since it does not take into account the 
energy dissipated (e.g. in shocks). Pcav is related to the 
pdV work through 

4PV 

Pcav — (2) 



3 



dChurazQv et al.|2QQ2[ pirzan et al.|2QQ8D , where r i 
duration of each sins;le AGN outburst and 4 is the f< 



is the 
factor 

used for relativistic plasma. Smolcic et al. ( 2009 ) investi- 
gated the fraction of radio AGM as a function of cosmic 
time and stehar mass of the galaxy. This fraction can 
be related, th r ough a probability argument detailed in 



Smolcic et al. (2009), with the time a galaxy of a given 
stellar mass and at a given redshift spends as a radio 
galaxy (rradio)- Using this result we can estimate the 
average duration of radio sources as a function of red- 
shift and stellar ma ss of the host galaxy (see Fig. 12 in 



Smolcic et al. 



2009) This gives a plausible time-scale 



during which the radio AGN can have injected mechan- 
ical energy into its environment. For the 16 COSMOS 
X-ray selected groups, Tradio ranges between 0.003 and 
4.18 Gyr, with a median value of 3.1 Gyr. The mechan- 
ical energy contribution can then be estimated as 



-E^mech — Pcav ^ '^radic 



(3) 



The values of Emech for our sources are shown in Ta- 
ble 1 and span a range between '^2x10^^-3x10^^ erg 
' The uncertainties in the radio mechanical energy 
input are dominated by the scatter in the scaling relation 
used to convert the monochromatic power into mechan- 
ical luminosity, which amounts to 0.85 dex, and by the 
uncertainties on Tradio- We use Tradio as derived from an 
average estimate over a sample of radio galaxies in the 
COSMOS field as a whole, irrespective of their environ- 
ment. One might expect the density of the environment 
surrounding the jets to have a significant impact on the 
jet lifetime. However, the fraction of radio galaxies that 
resides within the COSMOS groups is comparable with 
the fraction of red massive galaxies within groups in the 
control sample used in Smolcic et al. (2009) (respectively 
18% and 16% within R200); this assures that the statisti- 
cal argument used to compute the time-scales holds also 
in this case. Furthermore we can estimate an average 
time-scale based on only extended radio galaxies in the 
whole COSMOS group sample as follows. Of the 141 
COSMOS groups at ^ < 1 and with Lx > 10"^^ erg s"\ 
32 contain a multi-component radio galaxy. Therefore 
the average duration of the radio galaxy activity during 
this time interval is (32/141) x {t{z = 1) - t{z = 0)). 
This is ~ 1.7 Gyr, a time-scale comparable w ith the av- 
era ge life- time estimated with the method by [Smolcic et| 
^ (|2009[) (-1.6 Gyr^^ 



3.2. Binding energy of the intra-group medium 

We consider the shap e of the dar k matter halos to be 
characterized by NFW ( Navarro et al.„1996^ radial pro- 
files 



p{x) 



Pcrit 



(4) 



-•^^ To derive the values of Tradio I Smolcic et alj ( |2009| , it is as- 
sumed that the r adio par ent population (red massive galaxies) is 
formed at 2; = 3 ( |Renz ini 2006) and survives until z = 0. Since 
the COSMOS radio galaxies are not at 2; = 0, the time-scales com- 
puted in Smolcic et al. (2009) coincide with ours if multiplied by 
t(z—zgal)-t{z— 3.0) ^^lere t is the age of the universe at redshift z 

10^ yr ' ^ 

and zgal is the redshift of the radio galaxy. 



where x = r/rs^ Vg is the characteristic radius, and pc 
is the critical density of closure of the universe. 5c is 
defined as 

^ 200 
" 3 1n(l + c)-c/(l + c) 
and c is the concentration of the halo. The scale ra- 
dius and the concentration are linked by the relation 
Vs = i^5oo/c5oo, where C500 is the dark matter concen- 
tration inside i^soo- We estimate the binding energy out 
to i^5oo because the kinetic energy of the infall velocity 
field along filaments becomes important beyond this ra- 
dius (Evrard et al. 1996) and our simple model may not 
then be applicable. Furthermore, we can evaluate reli- 
able gas masses from the X-ray observations only within 
i?50o. For simplicity, we assume that the gas follows the 
same distribution as the dark matter. We define as bind- 
ing energy the total potential energy needed to push the 
ICM gas inside i?5oo beyond i^200- The binding energy 
is computed as 

^binding = / W) - ct){R200)] dMg 



--Air 



(j){r) pg{r) dr 



(6) 



We neglect the additive constant given by the term 
^(i?20o), as it is small with respect to the other terms of 
the equation. We use the definition of gas mass within 
R500 as 



M, = M,{R. 



500 j 



47r 



pg{r)r dr. (7) 



T he pot ential of a spherical NFW model is (Hayashi et 
^|2QQ7t 

, , , , ln(l + x) 
(I)(r) = Ax ^ (8) 



where A is 



200 



r,(/n(l + c)-c/(l + c)) 
Thus, substituting the terms into Equation |6| we com- 
pute the binding energy of the ICM gas in a NFW dark 
matter halo as follows: 

^binding = fgas^ TT Pcrit ^cAr^ I dx (10) 

Jo (J- + X) 

where /gas is the gas fraction. The concentration param- 
eter for the COSMOS groups has been computed from 



the mass-dependent relation of ^Maccio et al.| ( [2Q07[ ). The 



errors bars on the binding energy are estimated using a 
Monte Carlo method to numerically propagate the errors 
on M200 and i^200, the scatter in the C-M200 and in the 
/gas-^500 relation. 

We cannot estimate the gas masses from most of the 
existing X-ray observations of the COSMOS X-ray se- 
lected groups because of insufficient signal-to- noise. We 
therefore estimate the gas fraction in the groups from 
the mean trend of the gas mass fraction as a function of 
M500. This trend was established from an independent 
compilation of high quality observations of local (z<0.2) 
groups and clusters in the same ma ss range as th e sam- 



computed as 



max{T^(j^(ll(y) — min{r 



pie under consideration here (Pratt et al. 2009|). The 
observed relation (/gas ^'^tm) suggests that lower mass 
systems have proportionally less gas than high mass sys- 
tems. 
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TABLE 1 

For each of the 16 COSMOS groups, the columns indicate 1) X ray catalogue ID number 2) 
R.A. 3) Dec. 4) redshift 5) Power at 1.4 GHz 6) Mechanical Power 7)Binding energy ± la 

CONFIDENCE LIMIT 8) Emech ± CONFIDENCE LIMIT 9) Tradio 10) DISTANCE FROM THE CENTER 



XID 


R.A. 

[J2000] 


DEC. 
[J2000] 


z 


Pi AQ CtTJ 

[10^4 W/Hz] 


Pca-u 

[10^6 w] 


-^^b i n d i n g 

[lO^Oerg] 




[Gyr] 


R/R200 


107 


149.60965 


2.14799 


0.28 


1.11 


7.353 


p: qi ^2.342 


n c:-|43.087 
U.Oi^Q 44^ 


0.221 


0.3783 


262 


149.60007 


2.82118 


0.34 


19.0 


19.85 


-|27 Q38.80 
-^26.49 


94 1 cl44.9 
^^•-•-^20.70 


3.858 


0.0011 


253 


149.75626 


2.79472 


0.49 


6.71 


13.78 


40.16i2^if 


in 146O.85 


2.332 


0.0042 


246 


149.76132 


2.92909 


0.34 


0.90 


6.828 


289.8™:Ii 


7 nno42.02 


3.252 


0.6181 


311 


149.93796 


2.60627 


0.34 


6.38 


13.54 


17 974.792 

■■■'•^'3.907 


0-183j;i§7 


0.042 


0.2195 


264 


149.99847 


2.76914 


0.16 


0.32 


4.775 


6.940i:i™ 


n nn40.027 

U.UU^Q QQ3 


0.003 


0.0007 


281 


150.08617 


2.53141 


0.88 


8.90 


15.21 


84.94il;05 


1 9417.450 
-'-•^^-'-1.064 


0.258 


0.8614 


191 


150.11434 


2.35651 


0.22 


1.71 


8.554 


8.380i:?|9 


a 1 oo36.73 
"•-^^^5.248 


2.269 


0.0757 


237 


150.11774 


2.68425 


0.34 


27.7 


22.65 


io5.ii:!| 




2.738 


0.0027 


29 


150.17996 


1.76887 


0.34 


30.0 


23.29 


58.49i|:f0 




4.306 


0.0016 


64 


150.19829 


1.98628 


0.43 


12.4 


17.11 


91 in7.358 


10 n«l08.4 
iO.UO-L5 4g 


3.348 


0.0030 


35 


150.20661 


1.82327 


0.52 


10.2 


16.00 


on Q711.07 


1 c 0991.95 
-^^•'-*^13.13 


3.037 


0.0008 


6 


150.28821 


1.55571 


0.36 


1.13 


7.401 


77 Qi 20.23 
' ' -^-^17.35 


n n7n0.425 
^•^'^0.060 


0.030 


0.2279 


149 


150.41566 


2.43020 


0.12 


0.05 


2.564 


c:n 01 14.98 


n nnoO.014 

U.UUZq QQ2 


0.003 


0.1957 


40 


150.41386 


1.84759 


0.96 


48.5 


27.54 


iQo 134.85 
iUO.i24 40 


1 70^10.71 
'^"1.530 


0.205 


0.4888 


120 


150.50502 


2.22506 


0.83 


16.4 


18.88 


49c; nll2.5 


94 Q0I49.5 
^^•^^21.36 


4.185 


0.0267 



4. ANALYSIS OF GALAXY CLUSTER SAMPLE 

4.1. Mechanical energy input by radio galaxies in 
massive clusters 

In order to compare the energy input from radio galax- 
ies in groups and clusters, we include in our analysis a 
sample of well known radio gal axies in massive clus ters, 
extract ed fr om the sample of (Birzan et al. 2004 see 



Section 2.2). We use their tabulated value of pdV to 



compute the mechanical energy input over the average 
time the galaxy has spent as a radio galaxy. Birzan et 
al. provide a value for the energy input for both filled 
and radio ghost cavities. In order to obtain a measure 
of the average input, we sum the pdV for all the cavi- 
ties in a cluster and multiply it by the number of events 
(i.e. how often the radio jet was turned on). The latter 
is given by the ratio between Tradio and the duration of 
a single radio event (assuming that all the active AGN 
phases have the same duration). 

We choose the oldest cavity's ag e as an indica t ion of 
the duration of the radio event. Birzan et al. (2004) 
calculate the age of each cavity in three ways: 1^ tEe 
time required for the cavity to rise at the sound velocity; 
2. the time required for the bubble to rise buoyantly at 
the terminal velocity; 3. the time required to refill the 
displaced volume. We adopt the average of the three age 
estimates; this is generally similar to the age computed 
for a buoyantly rising bubble. We take the error on the 
cavity's age to be the difference between the shortest 
and longest life-time estimated via the three different 
methods. 

As T rad io foi" OUT Sample is derived following ISmolcic et 
al. ( 2009[ ), it depends on the redshift and the stellar mass 



of the radio galaxy. We computed the stellar masses for 
the central radio galaxy in the massive cluster sample us- 
ing the K-band photometry provided by the 2MASS sur- 
vey ( [Skrutskie et al.||2006[ ). This method is robust, since 
radio galaxies contam mostly type 2 (obscured) AGN, 
whose emission does not significantly contaminate the 
optical-NIR part of the galaxy spectrum. We assume a 
M/Lk ratio for a stellar population with an age of ^10 
Gyr (corresponding to the age of the stars in a galaxy at 



z-0), obtained by DroryetaL (2004') (M/Lk=1.4 with a 



Salpeter IM F). The quoted error on the M/Lk in Drory 
et al.| ( |2004 ) is 25-30%: a change in stellar mass of this 
magnitude does not affect significantly the time-scales 
we estimate. The stellar masses are then converted to a 
Chabrier IMF by subtracting an offset of 0.2 dex. 

4.2. Binding energy of the intra-cluster medium 

We compute the binding energy for the Birzan et al. 
clusters in the same way as for the COSMOS groups, 
using the value of M200 and R200 provided by the 
X-ray analysis in the HIFLUGCS survey (Reiprich & 
Boehringer 2002). We assume a constant concentration 
parameter of 5. Errors on £^binding are propagated nu- 
merically via a Monte Carlo method, in the same way 
as the COSMOS groups (see Section 3.2). As well as 
computing the binding energy of clusters individually, 
we also test the cluster result using the scaling relations 
adopte d for the COSMOS gr oups, both for computing 
M200 (Leauthaud et al.||201 0|) and for e stimating their 
mechanical energy output (Birzan et al.| |2008). The 
change in our calculations does not qualitatively affect 
our results. The values of £^binding change by less than a 



5 



factor 2 on average, while values of £^mech are perturbed 
randomly within the error bars. 

5. RESULTS 

5.1. The balance of radio-input and binding energy 

Figure [l] shows the binding energy of the gas versus 
the energy output from radio galaxies. In the group 
regime, the two energies span a comparable range of val- 
ues (10^^-10^^ ergs), while for clusters the binding energy 
exceeds the total mechanical output of radio-galaxies by 
a factor approximately of ~ 10^-10^. In particular, for 
seven groups the two energies are consistent at la level, 
and for all other groups except two the equality holds 
at 3(7, meaning that radio-galaxies potentially provide 
sufficient energy to unbind the gas in a large fraction 
of these groups. It is interesting to note that, in all the 
groups with ^mech ^^binding, the radio galaxy lies within 
0.15xR2oo from the center of the group. This suggests 
that a radio galaxy in a group is most likely to input 
sufficient energy into the ICM to unbind a part of the 
gas if it lies at the core of the group. Moreover, radio 
sources outside the group core reside in lower density en- 
vironments, and our calculations of those binding ener- 
gies may be overestimates. The different energy balance 
in groups and clusters demonstrates the importance of 
AGN heating in groups, and shows that the mechani- 
cal removal of gas from groups is energetically possible. 
This has important consequences for the understanding 
of the baryonic budget in these systems (see [Giodini et| 
aI] [2QQ9l ). 



5.2. Can radio galaxies offset radiative cooling in galaxy 

groups? 

We now compare the mechanical energy input by 
radio-galaxies with the energy required to offset the cool- 
ing in the g r oup center (^coo i)- As detailed in Fabian 



et al 
Nulsen 



(19941), IPeterson et al. d2003) and McNamara <H 



( [200^ 



nl^OT), the cooling time in cluster/group centers 



can be lower than the Hubble time, implying that large 
reservoirs of cold gas could accumulate in these regions. 
However, evidence that the gas does not cool b elow ap- 
proximately one third of the virial temperature (Kaastra 



et al. 2004) indicates the presence of a heat source pro- 
viding enough energy to offset the cooling. Several stud- 
ies (e.g. (Peterson et al. 2003 , IPeterson & Fabian||2006" 



[McNamara fc Nulsen||2007 ) suggest AGN feedback as a 
viable heating source. To test this hypothesis, we check 
whether the cooling energy is lower than the mechanical 
energy of the rising bubbles. We estimate ^coob assum- 
ing that the time during which the gas has been cooling 
is equal t o the lifetime of t he group, which we assume to 
be 5 Gyr ( Voigt fc Fabian|2004, ). The cooling energy can 
then be estimated as: 

-£^cool — -^cool ^ ty — /cool -^bol ^ ty {^^) 

where ty is the lifetime of the group, and /cool is the 
fraction of bolometric luminosity assumed to be emitted 
inside the cooling radius (where the cooling time of the 
gas is lower than the Hubble time). In general this con- 
tribution is fo und to be ^.10% of the to tal cluster X-ray 
luminosity (McNamara fc Nulsen |2007|). Also, the scat- 



ue to the contribution 



ter in the L^-T scaling relation c 
of cool core clusters can be up to a factor of 2 (iChen et 



LU 



10^' 



10^' 



~i 1 ;^T" 



O R<0.15R2oo 

O Clusters (Birzan et al. 2004) 
• COSMOS groups 



4- 



4 
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10"' 



Ebinding I^^Q] 



Fig. 1. — Binding energy of the cluster/group gas versus the 
output mechanical energy from radio-galaxies. Solid points show 
the 16 X-ray selected groups in the COSMOS field that host a radio 
galaxy within their virial radius. Large concentric circles indicate 
groups that host a radio galaxy within the core (R< 0.15 x i^20o)- 
Open points show the sample of massive local clusters drawn from 
[Birzan et al.|2004| The dashed line shows equality. The binding 
energy in clusters exceeds the total mechanical output by a factor 
of ~10-^-10^. In all cases except one where a radio galaxy lies in 
the centre of a group, the mechanical energy output from the radio 
galaxy is of the same order as the binding energy for the COSMOS 
groups analyzed here. 



al.||2QQ7 Pratt et al."2009). Given these considerations, 
we assume that 25% of the total bolometric X-ray lumi- 
nosity is emitted inside the cooling radius ( jPeres et al. 
1998). Since the relative contribution of the cool core 
to the total X-ray luminosity is higher in groups than 
in massive clusters, this value is a good estimate of the 
average contribution of the cooling core to the total lu- 
minosity of a group. 

In Figure [2] we compare ^mech and £^cooi in our groups. 
The mechanical energy injected by all but one of the core 
radio-galaxies is higher than the radiative losses, and ex- 
ceeds £^cooi by an order of magnitude in several cases. We 
can thus conclude that radiative losses do not greatly af- 
fect the net energy output of radio-galaxies in the cores 
of groups. On the other hand, the mechanical output by 
non-central radio galaxies is typically of the same order 
as £^cooi- Moreover, these sources reside mostly outside 
the cooling radius (~0.15 R2oo)5 where the cooling time 
is higher than the Hubble time. In this location, the gas 
does not lose as much energy through radiative cooling as 
in the core of the group, so these galaxies do not provide 
the required feedback at the right location. 

5.3. Impact of systematic effects 

The above calculations rest on several assumptions 
and should be regarded as rough estimates. One critical 
simplification is the calculation of the lifeti me of a 
radio-galaxy: the statistical argument used in |Smolcic 
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Fig. 2. — The output mechanical energy from radio-galaxies 
(Emech) versus energy radiated inside the cooling radius (Ecool] 
i.e. energy required to offset the cooling in the group center) for 
16 X-ray selected groups from the COSMOS survey . The dashed 
line shows the equality line. Red circles mark the radio-galaxies 
inside 0.15 x i^200- Uncertainties on Ecool cire computed allowing 
an error of a factor 2 on /cool- 



et al. (2009) relies on knowledge about the parent 
population that hosts the radio-galaxies. In the absence 
of evidence to the contrary, we assume that there is no 
significant difference between the radio-galaxy ellipti- 
cal hosts in groups and in l ow density environments 



(Feretti & Giovannini 2007). We note that even if 



'Tradio wcrc iucorrcct by a factor of 4, the mechanical 
output in clusters would still be significantly lower than 
the binding energy, but would remain consistent with 
the binding energy for many of the groups (see Figure [T]). 

Other biases may arise from the scaling relation of 
Birzan et al., which we use to compute the mechan- 
ical energy: the large scatter in the PiA9GHz~^cav 
relationship (0.85 dex) means that care must be taken 
when using the inferred value as the mean mechanical 
energy, since most of our calculations rest on the assump- 
tion that over the cluster/group lifeti me each burst has 
on ave rage the same power. Indeed, [Nipoti fc Binney] 
( 2005|) suggested that the distribution of the outbursts 
over the cluster /group lifetime is log-normal rather than 
gaussian; therefore in any system there would be a good 
chance of observing smaller than average jet powers. In- 
stead, much of the power would be generated by rare, 
more powerful outburst, s uch a s that observed in MS 
0735+7421 by |Gitti et aL] ( |2007D . These arguments rest 
on the assumption that the observed scatter in Pi. 49 ghz- 
Fcav in the observed ensemble of clusters is a good de- 
scription of the time variability of the AGN power in 
individual objects. In general the ensemble scatter is an 
upper limit to the scatter in the time variability. If we 
assume this scatter to represent also for the time variabil- 
ity, we are statistically underestimating the mechanical 
energy output over the group lifetime by a factor that we 
compute as follows. The scatter in the Birzan et al. re- 
lationship (0.85 dex) corresponds to a probability ^80% 
of observing a value smaller than the mean from a sin- 



gle observation (cf. [Nipoti fc Binney] 2005[ ). Thus, if 
we assume that the observed value of Piaoghz scatters 
around the median of the distribution, the ratio between 
the median and the mean for a lognormal distribution 
(which depends only on the scatter a) tells us the scal- 
ing factor for the 'true' mean mechanical energy: 



mean e 
median 2 



6.8. 



(12) 



Therefore the typical observed mechanical power may 
be underestimated by a factor <7 with respect to the 
mean. This value, though not negligible, goes in the di- 
rection of further increasing the mechanical output, con- 
firming the effect we found. 

Furthermore, if the bubble were ov er-pressured when 
compared to the surrounding ICM (Heinz et al. 1998), 
the expanding bubble would carry a shock and tne me- 
chanical p ower may be underestimated, as well as re- 
ported by [Birzan et al. (2004). This effect would also 
boost the mechanical energy to higher values, further 
strengthening our results. 

We have also used preliminary results from VLA 324 
MHz data (Smolcic et al. in preparation) to double- 
check our estimates of the mechanical energy output from 
radio-galaxies. Only 12 of the 16 radio-galaxies are de- 
tected in the 324 MHz band and, in all these case s, £^mech 
computed using these data (using Eq.l5 in Birzan et al. 
|2008[ ) is consistent within the error bars with the value 
computed at 1.49 GHz. As a further check, the total 
radio luminosity can be computed with higher precision 
fro m break frequencies for 7 of the 16 sources, using the 
the [Myers fc Spangler] ([1985|) approximation. The value 
of ^mech obtained with this improved method is con- 
sistent within the error bars with that obtained using 
monochromatic data. 

6. DISCUSSION: THE ENTROPY IN X-RAY GROUPS 

The injection of energy by radio galaxy activity into 
the ICM modifies the thermodynamical state of the gas, 
raising the entropy {S) by a significant amount compared 
to that generated by gravitational collapse. We define 
the entropy as 



kT 



(13) 



where T is the gas t emperature in keV and Ug is the 



gas electron density (Voit 



2005). An excess entropy 



of 50-100 keV cm~^ is indeed observed at the group 
regime, causing a deviation from the 5'-T relation 
(Ponman et al. 2003). The excess entropy is measured 
in the central regions (at 0.1 R2oo)- In the following we 
make an order of magnitude calculation of the excess 
entropy generated by the energy injected, then compare 
it with that observed in groups and predicted from the 
theory. We recompute the expected S-T relation taking 
into account this excess ene rgy, and compare it w ith the 
observational constraints of Ponman et al. (2003). 



The change in entropy caused by injection of energy 
under constant pressure is 



AS = 



2 AE 7^/^ 



5 ne 7t - 1 



(14) 



7 



( Lloyd-Davies et aL]|2QQQ ), where AE is the injected en- 
ergy per particle, 7t is the ratio between the init ial and 



final temperature (a value between 1.1 and 2.0, Lloyd 
Davies et al.||2Q QQ|), and ng is the initial electron density 
(ne=10~^ assuming the e nergy is deposited in the clus- 
ter core, e.g. ^Sanderson fc Ponman 2003) . We compute 
* " ' radio galaxies 

Ai? = i?„echX^ (15) 



from the mechanical energy input 
as follows: 



M, 



with M 



gas 



200' 



where 



'gas 



gas 

is estimated from the 



rela tion b etween gas fraction and total mass in |Pratt et 



al. (2009). This calculation does not depend on the de- 
tans oTthe energy injection process. We obtain values 
of excess entropy between 10 and 60 keV cm~^. This is 
a rough calculation but p redicts values similar to those 



Voit & Donahue (2005). These authors show that an 
additional energy input episodic on 10^ yrs timescale is 
needed to explain the excess entropy found observation- 
ally in the c ore of clusters ( [Ponman et al.|2003[ |Donahue| 
et al.|[2Q05|). The add itional energy produces an entropy 
pedestal: Voit (2005) calculates 10 keV cm~^ to be the 
minimum entropy boost needed to explain observations, 
and he predicts it to be larger for groups. 

The mechanical energy injected by radio galaxies 
into the 16 COSMOS X-ray selected groups is roughly 
independent on the group mass (see Figure [T]). This 
is not unexpected, since the black-hole masses (which 
are a ze roth order indicator of t he mechanical energy 
output; Merloni & Heinz 2007) range only between 
10^-10^ M(?) in radio galaxies (see Figure 7 in Smolcicet 



al. 2009). At the cluster regime may not be true that 
the me chanical energy is i ndependent on cluster mass. 
Indeed Chen et al. (2007) infer, from the strength of 
clusters' cooling cores, that a mechanical input higher 
than anything observed in groups is necessary to balance 
the cooling of the gas in the strong cool core clusters. 
However, it has been shown by the same authors that 
much (~90%) of that energy input would be radiated 
away to balance the cooling, and therefore would not 
participate to the mechanical removal of the gas. 

From these considerations, we can predict how the scal- 
ing relation between entropy and temperature is affected 
by the injection of a constant e xcess energy by radio 
galaxies. As shown in Finoguenov et all (120081), the en- 
ergy deposition into the ICM [AE) is proportional to 
the change in entropy ^ for a given typical Ug. We 
use ne=10~^ as the typical value of the density within 
0.1 R20O5 where th e majority of the energy is d eposited 



_2003). Using 

const ^ then 

(16) 



( deposition radius: Sanderson fc Po nman 
the scaling of Mgas ocT^ and Amech 

AS ^mech C 

where C is a constant and M^^s is the mass of the gas 
within the deposition radius. We can then infer the func- 
tional dependence of S on the virial temperature of the 
ICM as 



S = Sq^AS = Sqx{1 



AC n 

-) oc (To + (17) 



^0 ' 'To' 
where 5*0 and Tq are respectively the entropy and the 
temperature of the gas before the injection of energy from 



a radio galaxy . The value of C i s computed using Equa- 
tion 3 in Finoguenov et al. (2008 ) and has a median value 
of 2.56 if the energy is deposited inside the cooling ra- 
dius. We assunie the cooling radius to be 0.10 i?200 (e.g. 
Ponman et al.| '2003) and show the inferred functional 
form of t)(^E) in Figure [sj Remarkably, the shape of the 
resulting scaling relation (solid line) deviates from the 
self-similar one (dashed line) around ~4 keV, in agree- 
ment w ith the observed sca ling relation measured at 0.1 
1^200 by Ponman et al. ([2003 ) (black crosses; these points 
are binned means). The deviation of the '^1 keV point 
indicates that a lower excess entropy is needed to explain 
very cold groups. This can be achieved requiring that the 
mechanical energy is deposited at a larger radius in these 
groups. Indeed if the deposition radius increases also 
yigas within this radius increases. Therefore using Equa- 
tion 16 we would obtain a lower values of %^ (and thus 



entropy) for these groups, matching eventually the obser- 
vational point of Ponman et al. (2003) at ~1 keV; if this 
is the case, it would confirm that the eff ect of feedback is 
more global in groups than in clusters ( |De Young [2010 , 
cf.). Therefore, the injection of an excess energy tnat is 
independent of groups' mass, thus temperature, (as we 
observe from radio galaxies in the COSMOS groups) cor- 
rectly predicts the deviation of the observed S-T relation 
from the purely gravitational relation at the group scale. 



1000 




T, keV 



Fig. 3. — Scaling relation between entropy {S), measured at 0.1 
R200, and temperature (T). The solid black line is the inferred 
relation accounting for a constant energy excess injected by radio 
galaxies (see text for details). The grey line is the expected self- 
similar relation. The points show the binned means from the 
observations by (Ponman et al.| (|2003). The dashed line is the same 
as the solid line but considering the self-similar scaling M-gas ocT-*^-^ 



7. CONCLUSIONS 

In this Paper we have quantified the importance of the 
mechanical energy input by radio galaxies inside galaxy 
groups. In particular we report a striking difference be- 
tween clusters and groups of galaxies: while the binding 
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energy of the ICM in clusters exceeds the mechanical 
output by radio AGN, the two quantities are of the same 
order of magnitude in groups that host a radio galaxy 
within 0.15 R2oo- This suggests that, while clusters can 
be mostly considered to be closed systems, the mechan- 
ical removal of gas is energetically possible from groups. 
This has implications that help explain recent findings 
on the bary onic fraction in groups of galaxies. Giodini 



et al. 



( Dunkley et al.|[2QQ9| . It has been suggested that 
k has been removed by AGN feedback. 



( |2QQ9[ ) reported a ~30% lack of gas m groups com- 
pared with the cosmological baryon mass fraction evalu- 
ated fr om the 5 years Wilk inson Microwave Anisotropy 
Probe 

this gas iias been removed by 

This is consistent with cosmological models in which 
feedback from radio galaxies is invoked to successfully 
explain galaxy group/cluster properties. Based on a well 
selected sample of galaxy groups and clusters that host 
radio galaxies, we have observationally shown for the 
first time that this scenario is energetically feasible. We 
have further shown that a constant injection of excess en- 



ergy by radio galaxy naturally reproduces the self-similar 
breaking observed in the scaling relation between the en- 
tropy and temperature of groups. 
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APPENDIX 
RADIO IMAGES 



Fig. Al. — Figure [aT] presents the contours maps of the radio 20 cm emission (magenta hues) superimposed to the SUBARU zp band 
image for each of the groups hsted in Table Images are 3x3 arcmin and centered on the group center, except XID246 which is 4x4 
arcmin wide and offset from the center group oecause of its location on the edge of the SUBARU field coverage. The white line shows the 
contours of X-ray fiux significance. The contours correspond to [3, 6, 9, 12, 15, 18, 21, 24] a X-ray fiux significance. 
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